Parametric autoresonance 
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We investigate parametric autoresonance: a persisting phase locking which occurs when the driving 
frequency of a parametrically excited nonlinear oscillator slowly varies with time. In this regime, the 
resonant excitation is continuous and unarrested by the oscillator nonlinearity. The system has three 
characteristic time scales, the fastest one corresponding to the natural frequency of the oscillator. We 
perform averaging over the fastest time scale and analyze the reduced set of equations analytically 
and numerically. Analytical results are obtained by exploiting the scale separation between the two 
remaining time scales which enables one to use the adiabatic invariant of the perturbed nonlinear 
motion. 

PACS numbers: 05.45.-a 



I. INTRODUCTION 

This work addresses a combined action of two mech- 
anisms of resonant excitation of (classical) nonlinear os- 
cillating systems. The first is parametric resonance. The 
second is autoresonance. 

There are numerous oscillatory systems which interac- 
tion with the external world amounts only to a periodic 
time dependence of their parameters. The corresponding 
resonance is called parametric [QJ^] . A textbook example 
is a simple pendulum with a vertically oscillating point of 
suspension jlj] . The main resonance occurs when the ex- 
citation frequency uj is nearly twice the natural frequency 
of the oscillator uiq [0j2| . Applications of this basic phe- 
nomenon in physics and technology are ubiquitous. 

Autoresonance occurs in nonlinear oscillators driven by 
a small external force, almost periodic in time. If the 
force is exactly periodic, and in resonance with the nat- 
ural frequency of the oscillator, the resonance region of 
the phase plane has a finite (and relatively small) width 
HI). If instead the driving frequency is slowly varying 
in time (in the right direction determined by the nonlin- 
earity sign), the oscillator can stay phase- locked despite 
the nonlinearity. This leads to a continuous resonant ex- 
citation. Autoresonance has found many applications. It 
was extensively studied in the context of relativistic par- 
ticle acceleration: in the 40-ies by McMillan ||, Vcksler 
and Bohm and Foldy [fTp], and more recently [p|-p2[. 
Additional applications include a quasiclassical scheme of 
excitation of atoms |1S j and molecules [14j], excitation of 
nonlinear waves jl5|, 16 1, solitons |l7|]l§| ], vortices [ p~9|j20) 
and other collective modes |2l|] in fluids and plasmas, an 
autoresonant mechanism of transition to chaos in Hamil- 
tonian systems [p2[^3| , etc. 

Until now autoresonance was considered only in sys- 
tems executing externally driven oscillations. In this 
work we investigate autoresonance in a parametrically 
driven oscillator. 

Our presentation will be as follows. In Section 2 we 
briefly review the parametric resonance in non-linear os- 



cillating systems. Section 3 deals, analytically and nu- 
merically, with parametric autoresonance. The conclu- 
sions are presented in Section 4. Some details of deriva- 
tion are given in Appendices A and B. 



II. PARAMETRIC RESONANCE WITH A 
CONSTANT DRIVING FREQUENCY 

The parametric resonance in a weakly nonlinear oscil- 
lator with finite dissipation and detuning is describable 
by the following equation of motion [p 24 25f|: 



x + 2-^x + [1 + e cos ((2 + 6)t)} x - fix 3 



0. 



(1) 



where the units of time are chosen in such a way that the 
scaled natural frequency of the oscillator in the small- 
amplitude limit is equal to 1. In Eq. (|]) e is the ampli- 
tude of the driving force, which is assumed to be small: 
0<e<1, (S<lis the detuning parameter, 7 is 
the (scaled) damping coefficient (0 < 7 < 1) and (3 is 
the nonlinearity coefficient. For concreteness we assume 
(3 > (for a pendulum (3 = 1/6). 

Working in the limit of weak nonlinearity, dissipa- 
tion and driving, we can employ the method of aver- 
aging P,p p^j27| ], valid for most of the initial conditions 
HQ- The unperturbed oscillation period is the fast time. 
Putting x = a(t) cos Q(t) and x = —a(t) sin#(i) and per- 
forming averaging over the fast time, we arrive at the 
averaged equations 
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where a new phase if> = — [(2 + S)/2]t has been in- 
troduced. The averaged system (^) is an autonomous 
dynamical system with two degree of freedom and there- 
fore integrable. In the conservative case 7 = Eqs. (Q) 
become: 
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(3) 



As sin 2ip and cos 2ip are periodic functions of ip with a 
period 7r, it is sufficient to consider the interval —ir/2 < 
ip < tt/2. For small enough detuning, 5 < e/2, there is 
an elliptic fixed point with a non-zero amplitude: 
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ip* = o. 



We need to calculate the period of motion in the phase 
plane along a closed orbit around this fixed point (such 
an orbit is shown in Fig. [I]). 
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FIG. 1. Parametric resonance with a constant driving fre- 
quency. Shown is a typical closed orbit in the phase plane 
with period T (e = 0.04, S = -0.04, (3 = 1/6 and 7 = 0.). For 
a time-dependent driving frequency v{t), the autoresonance 
will occur if the characteristic time for variation of v(t) is 
much greater than T, see criterion (|l"l|). 

This calculation was performed by Struble 0|. For a 
zero detuning, 5 — 0, Hamilton's function (we will call it 
the Hamiltonian) of the system (|J) is the following: 

el 30 1 2 

H(J, ip) = -r cos 2ip - —— =H = const., (4) 
4 8 

where we have introduced the action variable / = a 2 /2. 
Solving Eq. (Q) for / and substituting the result into the 
Hamilton's equation for ip we obtain: 
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(5) 



where the minus (plus) sign corresponds to the upper 
(lower) part of the closed orbit. The period of the am- 
plitude and phase oscillations is therefore 



T = 



dip 



(cos 2 2ip 
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where — ip and ip are the roots of the equation cos 2 2ip 
24/3 Ho/e 2 . Calculating the integral, we obtain 



T 



-K(m) 



(7) 



where K(m) is the complete elliptic integral of the first 
kind @, and m = 1 - 2A0H o /e 2 . This result will be 
used in Section 3 to establish a necessary condition for 
the parametric autoresonance to occur. 



III. PARAMETRIC RESONANCE WITH A 
TIME-DEPENDENT DRIVING FREQUENCY: 
PARAMETRIC AUTORESONANCE 

Now let the driving frequency vary with time. This 
time dependence introduces an additional (third) time 
scale into the problem. The governing equation becomes 



+ 2 7 i; + (1 + 
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j3x 3 = 0, 



(8) 



where cp = v(t). We will assume v{t) to be a slowly de- 
creasing function which initial value is v(t = 0) = 2 + 5. 
Using the scale separation, we obtain the averaged equa- 
tions. The averaging procedure of Section 2 can be re- 
peated by replacing (2 + S)t by (p in all equations. There 
is one new point that should be treated more accurately. 
The averaging procedure is applicable (again, for most 
of the initial conditions) if there is a separation of time 
scales. It requires, in particular, a strong inequality 
20 + v(t) ^ 26 — v(t). This inequality can limit the time 
of validity of the method of averaging. Let us assume, 
for concreteness, a linear frequency "chirp" : 



v(t) = 2 + 5 - 2/jt, 



(9) 



where fj, <C 1 is the chirp rate. In this case the averaging 
procedure is valid as long as [it ^ 1. 

Introducing a new phase ip = 9 — 4>/2, we obtain a 
reduced set of equations (compare to Eqs. Q): 
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(10) 



The first of Eqs. ( |Io|) is typical for parametric resonance: 
to get excitation one should start from a non-zero os- 
cillation amplitude. As we will see, the fit term in the 
second of Eqs. (|lTj) (when small enough and of the right 
sign) provides a continuous phase locking, similar to the 
externally driven autoresonance. 

Consider a numerical example. Fig. E shows the time 
dependence a(t) found by solving Eqs. (10) numerically. 
One can see that the system remains phase locked which 
allows the amplitude of oscillations to increase, on the 
average, with time in spite of the nonlinearity. The time- 
dependence of the amplitude includes a slow trend and 
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relatively fast, decaying oscillations. These are the two 
time scales remaining after the averaging over the fastest 
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FIG. 2. An example of parametric autoresonance. Shown 
is the oscillation amplitude versus time, computed nu- 
merically from the averaged equations (hoj) . The system 
remains phase-locked which allows the amplitude to in- 
crease, on the average, with time. The parameters are 
H = 6.5 ■ 1(T 5 , e = 0.04, 5 = -0.01, (3 = 1/6 and 7 = 0.001. 

Similar to the externally-driven autoresonance, a per- 
sistent growth of the oscillation amplitude requires the 
characteristic time of variation of v(t) to be much greater 
than the "nonlinear" period T [see Eq. (Q)] of oscillations 
of the amplitude: 



(11) 



Like its externally-driven analog, the parametric au- 
toresonance is insensitive to the exact form of v(t). For 
a given set of parameters, the optimal chirping rate can 
be found: too low a chirping rate means an inefficient ex- 
citation, while too high a rate leads to phase unlocking 
and termination of the excitation. 

In the remainder of the paper we will develop an an- 
alytical theory of the parametric autoresonance. The 
first objective of this theory is a description of the slow 
trend in the amplitude (and phase) dynamics. When the 
driving frequency v is constant, there is an elliptic fixed 
point a* (see Section 2). When v varies with time, the 
fixed point ceases to exist. However, for a slowly-varying 
v(t) one can define a "quasi-fixed" point a*(t) which is 
a slowly varying function of time. It is this quasi-fixed 
point that represents the slow trend seen in Fig. |^ and 
corresponds to an "ideal" phase-locking regime. The fast, 
decaying oscillations seen in Fig. ^ correspond to oscil- 
lations around the quasi-fixed point in the phase plane 
[this phase plane is actually projection of the extended 
phase space (a,ip,t) on the (a, ip) -plane]. 

In the main part of this Section we neglect the dissi- 
pation and use a Hamiltonian formalism. First we will 
consider excitation in the vicinity of the quasi-fixed point. 



Then excitation from arbitrary initial conditions will be 
investigated. Finally, the role of dissipation will be briefly 
analyzed. 

For a time-dependent v(t), the Hamiltonian becomes 
[compare to Eq. (0)] : 



H(I,iP,t) = j (a(t) + cos2V>) - 



(12) 



where a(t) = (4/e)(l — z/(i)/2). The Hamilton's equations 
are: 
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Let us find the quasi-fixed point of (|l3|), i.e. the spe- 
cial autoresonance trajectory I*(t), ip*{i) corresponding 
to the "ideal" phase locking (a pure trend without oscil- 
lations). 

Assuming a slow time dependence, we put ip* — 0, that 

is 



— (a + cos 2ip* ) ; — = 0. 

4 4 



(14) 



Differentiating it with respect to time and using Eqs. 
(O), we obtain an algebraic equation for ip„(t): 



2a(t) sin 2^* + sin 4^* 



16(i 



(15) 



At this point we should demand that ip* (t) , evaluated on 
the solution of Eq. (|l5|), is indeed negligible compared 
to the rest of terms in the equation dl3| ) for ijj{t). It is 
easy to see that this requires 16/^/e 2 < 1. In this case 
the sines in Eq. (|l5| ) can be replaced by their arguments, 
and we obtain the following simple expressions for the 
quasi-fixed point: 



(16) 



where k = 4///e 2 



A. Excitation in the vicinity of the quasi-fixed point 

Let us make the canonical transformation from vari- 
ables / and i/j to SI = I — /* and Sip — ip — ip* . Assuming 
81 and Sip to be small and keeping terms up to the second 
order in SI and Sip, we obtain the new Hamiltonian: 



H{SI,Sip,a(t)) = 



a + l 



SISip 



3 



Here and in the following small terms of order of k 2 are 
neglected. Let us start with the calculation of the local 
maxima of SI(t) and 5ip(t), which will be called 6I max (t) 
and Sip max (t), respectively. As a(t) is a slow function of 
time [so that the strong inequality ( pd| ) is satisfied], we 
can exploit the approximate constancy of the adiabatic 
invariant [02911: 



J = 



1 

2^ 



SId{5tp) ~ const. 



(18) 



| J | is the area of the ellipse defined by Eq. ([17]) with the 
time-dependencies "frozen". Therefore, 
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const. 



(19) 



e (a+1) 1 / 2 

This expression can be rewritten in terms of SI and 5tp: 
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+ ^(a + l) 1/ W- 



(20) 



If fc = 4/i/e 2 <C 1, the term with SISif) in (gq) can be 
neglected (in this approximation one has ip* = 0). Then 
J becomes a sum of two non- negative terms, one of them 
having the maximum value when the other one vanishes. 
Therefore, 



SIraax (t) 2 1 



and 



84>max(t) 



= ( 3/3J ) 



1 
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(22) 



Now we calculate the period of oscillations of the ac- 
tion and phase. Using the well-known relation fjj T = 
2-K{dJ/dH), we obtain from Eq. @: 
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The period of oscillations versus time is shown in Fig. 
The theoretical curve [Eq. (23)] shows an excellent agree- 
ment with the numerical solution. 

Now we obtain the complete solution 8I(t) and 5ip(t). 
The Hamilton's equations corresponding to the Hamilto- 
nian (1l7|) are: 



SI 

Sip : 
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(24) 



Differentiating the second equation with respect to 
time and substituting the first one, we obtain a linear 
differential equation for Si/j(t): 



S'^ + uj 2 {t)S^ = 0, 



(25) 



where u(t) = (e/2)(a(t) + 1) 1/2 . For the linear v(t) de- 
pendence (Eq. (ph) we have a(t) = Aut/e — 25/e, therefore 
for C 1 the criterion lu/lo 2 <gc 1 is satisfied, and Eq. 
( p5| ) can be solved by the WKB method (see, e.g. [[|). 
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FIG. 3. Excitation in the vicinity of the quasi-fixed point: 
the time-dependence of the period T of the action and phase 
oscillations. The solid line is the theoretical curve, Eq. (p3[), 
the asterisks are points obtained numerically. The parameters 
are fi = 6.5 ■ 1(T 5 , e = 0.04, 5 = -0.01 and /3 = 1/6. 

The WKB solution takes the form (details are given in 
Appendix A): 



V e J (a + l)V4 
(a + 1) 3 / 2 
3fc 



X COS q + 



(26) 



where the phase qo is determined by the initial conditions. 
The full solution for the phase is ip = Sip + ip* and Fig. ^ 
compares it with a numerical solution of Eq. jl^). Also 
shown are the minimum and maximum phase deviations 
predicted by Eqs. ( p2|) and (|l^). One can see that the 
agreement is excellent. 

The solution for 51{t) can be obtained by substituting 
Eq. ( pif ) into the second equation of the system (^4|). I n 
the same order of accuracy (see Appendix A) 



«(i)=2^y /2 (a + l) 1/4 sin L + 



(a + 1) 3 / 2 



3A- 



(27) 



Fig. |^ shows the dependence of the action variable 
with the trend I*(t) subtracted, 8I(t), on time predicted 
by Eq. (p7|), and found from the numerical solution. It 
also shows the minimum and maximum action deviations 
(pl|). Again, a very good agreement is obtained. 
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FIG. 4. Parametric autoresonance excitation in the vicin- 
ity of the quasi-fixed point. Shown is the phase ip(t) found 
analytically [Eqs. (|l^) and ([S])] and by solving Eq. ( |is| ) nu- 
merically. The analytical and numerical curves are indistin- 
guishable. Also shown are the minimum and maximum phase 
deviations predicted by Eq. (^ij) and (|l6|). The parameters 
are the same as in Fig. ||| 
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FIG. 5. Parametric autoresonance excitation in the vicinity 
of the quasi-fixed point. Shown is the action variable SI(t) 
from Eq. ( p7| ) and from the numerical solution. Also shown 
are the minimum and maximum action deviations predicted 
by Eq. (^). The parameters are the same as in Fig. ^. 



B. Excitation from arbitrary initial conditions 

In this Subsection we go beyond the close vicinity oi 
the quasi-fixed point and calculate the maximum devi- 
ations of the action / and phase ^ for arbitrary initial 
conditions. Again, these calculations are made possible 
by employing the adiabatic invariant for the general case. 
Correspondingly the period of the action and phase os- 
cillations will be also calculated. 

Let us first express the maximum and minimum actior, 
deviations in terms of the Hamiltonian H and driving 
frequency v{t). Solving Eq. (|l2|) as a quadratic equation 



for /, we obtain: 



h.2 = 3^ ( a + cos2 V0 ± 



e 2 , „ ,s2 8# 

9^(a + cos2V0 -30 



1/2 



The time derivative of / vanishes when / = I max or 
I = Imin- Therefore, from the first equation of the sys- 
tem (O) V = so that 



' HI (t.r.llt i li 



3/3 



(a + l)± 
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9/3~ 



£ / , \2 $>H U p down 

2 (a + 1) 



3f3 



1/2 



(28) 



where H U p a ' owrL H(Imax,min>ip — 0). 

Now we express the maximum and minimum phase 
deviations through the Hamiltonian H and driving fre- 
quency v{t). The time derivative ip vanishes if ip = ipmax 

Or I p — Ipmin 5 

then the second equation of the system 
( |l3| ) yields / = (e/3/3)(a + cos 2-0). In this case the 
Hamiltonian (|l2|) becomes H rig i lt ,ieft = (e 2 /24/?)(a + 



cos 2tp. 



max , min 



) 2 . Finally, the expression for ip. 
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IS 
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±— arccos 
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(29) 



Fig. H shows a part of a typical autoresonant orbit 
in the phase plane. For v{t) — const, this orbit is de- 
termined by the equation H^I^ip^v) — const. , and it is 
closed. As in our case v{t) changes with time, the tra- 
jectory is not closed. To calculate the maximum and 
minimum deviations of action and phase we should know 
the values of the Hamiltonian at 4 points of the orbit 
that we will call "up" , "down" , "left" , and "right" in the 
following. 
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FIG. 6. A part of the autoresonant orbit in the phase plane. 
Knowing the Hamiltonian at the 4 points, we can calculate the 
maximum and minimum deviations of the action and phase. 
The parameters are the same as in Fig. ^. 
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Knowing the values of the Hamiltonian at these 4 
points, we calculate I m ax,min horn Eq. (|2f) and ip m ax,min 
from Eq. (|29|). Figs. (0) and (||) show these deviations 
for action and phase correspondingly, and the values of 
/ and ijj, found from numerical solution. The theoretical 
and numerical results show an excellent agreement. 
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FIG. 7. The maximum and minimum deviations of the ac- 
tion, calculated from Eq. (thick line) and from numerical 
solution (thin line). The parameters are the same as in Fig. 
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FIG. 8. The maximum and minimum deviations of the 
phase, calculated from Eq. (^) (thick line) and from numer- 
ical solution (thin line). The parameters are the same as in 
Fig. | 



so that Eq. (Bffl) becomes: 



J = — / arccos a ) dl, (32 



where I max and I m in are given by Eq. (|2q). Notice that 
H(t) and a(t) should be treated as constants under the 
integral (|3|), see Refs. [|l]]|,|2j|. This integral can be ex- 
pressed in terms of elliptic integrals (see Appendix B for 
details). For definiteness, we used the values of H{t) 
and a(t) in the "up" points, see Fig. |(| We checked 
numerically that the adiabatic invariant J(H(t),a(t)) is 
constant in our example within 0.12 per cent. 

Now we calculate the period of action and phase os- 
cillations. From the first equation of system (|l^) we 
have: 



T = 2 



dl 



(e//2)sin2V>' 



(33) 



where I max and I mm are given by Eq. @, while ib 
ib(I) is defined by @. 

Using Eq. (p"2|), we obtain after some algebra: 



T = — 



dl 



30J Imin G(I) 1 / 2 ' 



(34) 



where G(I) is given in Appendix B, Eq. (B2). Again, we 
treat H (t) and a(t) as constants under the integral (|34|) , 
and take their values in the "right" points, see Fig. |6|. 
The final result is: 



where Co 



T = C7 2 A-(C7 3 ), 
4(2/3/3i/e 2 ) 1 / 4 and 



(35) 



C 3 



CI 
16 



3/3H 



16 



(1 



Now we are prepared to calculate the adiabatic invari- 
ant J(H,v(t)). Its (approximate) constancy in time al- 
lows one, in principle, to find the Hamiltonian H(t) at 
any time t, in particular at the points of the maximum 
and minimum action and phase deviations (see Fig. ^J). 

It is convenient to rewrite the adiabatic invariant in 
the following form: 



Using Eq. @, we can find ip = ip{H, I, a(t)): 
'8H + 3/3I 



ib = ±— arccos 
r 2 



2el 



(30) 



(31) 
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FIG. 9. The period T of the phase (action) oscillations ob- 
tained from Eq. ( pS5| ) (solid line), and from numerical solution 
(asterisks). The parameters are the same as in Fig. |j| 
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Figure ^ shows the period T of the phase and action 
oscillations versus time obtained analytically and from 
numerical solution. This completes our consideration of 
the parametric autoresonance without dissipation. 



C. Role of dissipation 

Now we very briefly consider the role of dissipation 
in the parametric autoresonance. Consider the averaged 
equations ([lO]) and assume that the detuning is zero. The 
non-trivial quasi-fixed point exists when the dissipation 
is not too strong: 7 < e/4, and it is given by 



2e 



1/2 



ip* = — arcsin 



a(t) 
47 



f - 



1/2 



2k 



a(i) + (f-f6 7 2 /e 2 ) 1 /2 



(36) 



Again, we assume k <C 1. This quasi- fixed point de- 
scribes the slow trend in the dissipative case. As we 
see numerically, fast oscillations around the trend, 5a — 
a — a* and Sip — ip — ip* decay with time. Therefore, one 
can expect that the a(t) will approach, at sufficiently 
large times, the trend a*(t). Fig. |l^ shows the time de- 
pendence of the amplitude, found by solving numerically 
the system of averaged equations (|lCj) , and the amplitude 
trend from (|36|). We can see that indeed the amplitude 
a(t) approaches the trend a^(t) at large times. 
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FIG. 10. Parametric autoresonance with dissipation: the 
time dependence of the amplitude of oscillations, obtained 

from numerical solution of Eqs. ( lOj) , and the amplitude 

trend a*(t), predicted by Eq. (p6). The parameters are 
M = 6.5 • 1(T 5 , e = 0.04, 8 = 0, 7 = 0.002 and (3 = 1/6. 



IV. CONCLUSIONS 

We have investigated, analytically and numerically, a 
combined action of two mechanisms of resonant excita- 
tion of nonlinear oscillating systems: parametric reso- 
nance and autoresonance. We have shown that para- 
metric autoresonance represents a robust and efficient 
method of excitation of nonlinear oscillating systems. 
Parametric autoresonance can be extended for the ex- 
citation of nonlinear waves. We expect that parametric 
autoresonance will find applications in different fields of 
physics. 
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APPENDIX A: CALCULATION OF PHASE AND 
ACTION DEVIATIONS BY THE WKB-METHOD 

Changing the variables from time t to a, we can rewrite 
Eq. (p5| ) in the following form: 



„ , a (t) + 1 
4k 2 



Sip" + 



(Al) 



where " denotes the second derivative with respect to a. 
Solving this equation by the WKB-method j| , we obtain 
for dtp: 

wt) .^ M fa + («('>+.y- 1 V (A2) 



(a + l) 1 ^ — r" u ' 3fc 
where £Iq and C are constants to be found later. Now 



we obtain the solution for SI. Substituting (A2) into the 
second equation of the system (^4|) , we obtain in the same 
order of accuracy: 



SI(t)^^(2kC)^(a + l)^ 



x sin fl 



in 



(a(t) + If' 2 - 1 



3k 



(A3) 



The constant C can be expressed thro ugh the adiabatic 
invariant J, given by (|20|). From Eqs. (A2) and (A3) we 
have: 



Therefore, a small amount of dissipation enhances 
the stability of the parametric autoresonance excitation 
scheme. A similar result for the externally-driven au- 
toresonance was previously known |p0[. 



Comparing it with ( pTj| ) we find: C ~ 3/3J/2fce. Substi- 
tuting this value into Eqs. (|A2|) and ( jA3[ ) we obtain the 
final expressions (pq) and (|27|) for Sip(t) and SI(t). 
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APPENDIX B: CALCULATION OF THE 
ADIABATIC INVARIANT 

After integration by parts and some algebra, using 
Eqs. ([[2]) and (|28|), we obtain the following expression 
for the adiabatic invariant: 



J 



1 

2^ 



' T 2 SH ' 
1 3/3 



Imin 



Gil) 1 / 2 



(II. 



where 

G(I) = (/ maa -/)(/- 



e(l-Q) 
3/3 



(Bl) 



16D 



9/3 2 
(B2) 



and we assume D = (e 2 /16)(l - a) 2 - 3/3ff/2 < 0. Calcu- 
lation of this integral employs several changes of variable 
shown in the best way by Fikhtengolts piR . Using the 
reduction formulas [B8|, we arrive at: 



J = C X 
1 



1 + mm! 



n 



(1 — m) 2 (l + to') \m — 1 



m 



-\k' 



1 — TO 



K(k 2 ) 



to + m! 



where 



(1 - m)(l + to') 
(e/3/3)(l + a) + (8 J ff/3/3) 1 /2 



(B3) 



>0, 



, (elZm-a) + m/Zf3) 1/2 ^ n 

777 =^ ~> II 

-(e/3/3)(l-a) + (8H/3py/ 2 



TO + TO' 



3/3(m + m') 1/2 ' 



and 



1 e /8/f 
2^ 3^ (1 + a)+ U 



1/2 



-1/2 



'3^ (1 - a)+ U 
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-1/2 



[2 



[3. 
K 
[5 

ir 

[8 
[9 
[10 

[11 

[12. 
[13 

[14 
[15 

[16 
[17 

[is; 

[19 
[20 

[21 

[22 
[23 
[24 
[25 : 
[26 



[27 

[28 

[29 

[30 
[31 



N.N. Bogolubov and Y.A. Mitropolsky, Asymptotic 
Methods In The Theory of Non-linear Oscillations (Gor- 
don and Breach Science Publishers, New York, 1961). 
R.Z. Sagdeev, D.A. Usikov, and G.M. Zaslavsky, Nonlin- 
ear Physics (Harwood Academic, Switzerland, 1988). 

A. J. Lichtenberg and M.A. Lieberman, Regular and 
Chaotic Dynamics (Springer- Verlag, Oxford, 1992). 
E.M. McMillan, Phys. Rev. 68, 143 (1945). 

V. Veksler, J.Phys.(USSR) 9, 153 (1945). 
D. Bohm and L. Foldy, Phys. Rev. 70, 249 (1947). 
D. Bohm and L. Foldy, Phys. Rev. 72, 649 (1947). 
K.S. Golovanivsky, Phys. Scripta 22, 126 (1980). 

B. Meerson, Phys. Lett. A 150, 290 (1990). 

B. Meerson and T. Tajima, Optics Communications 86, 
283 (1991). 

L. Friedland, Phys. Plasmas 1, 421 (1994). 

B. Meerson and L. Friedland, Phys. Rev. A 41, 5233 

(1990). 

J.M. Yuan and W.K. Liu, Phys. Rev. A 57, 1992 (1998). 
M. Deutsch, B. Meerson, and J.E. Golub, Phys. Fluids 
B 3, 1773 (1991). 

L. Friedland, Phys. Plasmas 5, 645 (1998). 

I. Aranson, B. Meerson, and T. Tajima, Phys. Rev. A 

45, 7500 (1992). 

L. Friedland and A. Shagalov, Phys. Rev. Lett. 81, 4357 
(1998). 

L. Friedland, Phys. Rev. E 59, 4106 (1999). 

L. Friedland and A.G. Shagalov, Phys. Rev. Lett. 85, 

2941 (2000). 

J. Fajans, E. Gilson, and L. Friedland, Phys. Rev. Lett. 

82, 4444 (1999); Phys. Plasmas 6, 4497 (1999). 

B. Meerson and S. Yariv, Phys. Rev. A 44, 3570 (1991). 

G. Cohen and B. Meerson, Phys. Rev. E 47, 967 (1993). 
R.A. Struble, Quart. Appl. Math. 21, 121 (1963). 
A.D. Morozov, J. Appl. Math. Mech. 59, 563 (1995). 
M.I. Rabinovich and D.I. Trubetskov, Oscillations and 
Waves in Linear and Nonlinear Systems (Kluwer Aca- 
demic Publisher, Dordrecht, 1989). 

P.G. Drazin, Nonlinear Systems (Cambridge University 
Press, Cambridge, 1992). 

M. Abramowitz, Handbook of Mathematical Functions 
(National Bureau of Standards, Washington, 1964). 

H. Goldstein, Classical Mechanics (Addison- Wesley, 
Reading, Mass., 1980). 

S. Yariv and L. Friedland, Phys. Rev. E 48, 3072 (1993). 
G.M. Fikhtengolts, The Fundamentals of Mathematical 
Analysis (Pergamon Press, New York, 1965). 



Here K, E and n are the complete elliptic integrals of 
the first, second and third kind, respectively. 



[1] L.D. Landau and E.M. Lifshits, Mechanics (Pergamon 
Press, Oxford, 1976). 



8 



